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We show the role played by QCD evolution and skewedness effects in the DVCS cross section at 
large Q 2 within the color dipole description of the process at photon level. The dipole cross section 
is given by the saturation model, which can be improved by DGLAP evolution at high photon 
virtualities. We investigate both possibilities as well as the off-forward effect through a simple 
phenomenological parametrisation. The results are compared to the recent ZEUS DVCS data. 
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I. INTRODUCTION 

An important clean process allowing us to access off- 
diagonal (skewed) parton distributions, which carry new 
information on the nucleon's dynamical degrees of free- 
dom, is the Deeply Virtual Compton Scattering (DVCS) 
IH 13 This is due to the real photon in the final state 
being an elementary (point-like) particle rather than a 
bound state like a meson or more complicated configu- 
rations. The skewed parton distributions are generally 
defined via the Fourier transform of matrix elements of 
renormalized, non-local twist-two operators (for a peda- 
gogical view, see Refs. Hi). These composite operators 
contain only two elementary fields of the theory, which 
are placed at different positions becoming then non-local 
and operating in unequal momentum nucleon states. 

Hence skewedness takes into account dynamical cor- 
relations between partons with different momenta. The 
high energy situation at HERA gives the important op- 
portunity to constrain them as well as to study the evolu- 
tion with virtuality of the resulting quark and gluon dis- 
tributions. There are several representations for skewed 
parton distributions [ESSE! EH > which can be used 
to compute the relevant observables in DVCS (or other 
exclusive processes) through a factorization theorem . 
They are input in numerical solutions of the renormal- 
ization group or evolution equations (see e.g. |l2 l). p ro- 
ducing very reliable predictions up to NLO level [Lj 

On the other hand, the color dipole models have 
also been successful in describing DVCS observables 
(Til ITU Ht| . There, the main degrees of freedom are 
the color dipoles, which interact with the nucleon tar- 
get via gluonic exchange. This interaction is modeled 
through the dipole-nucleon cross section, which can in- 
clude QCD dynamical effects given by DGLAP, BFKL or 
non-linear high energy evolution equations (parton sat- 



uration). Skewedness effects are not considered in the 
current dipole models and this is one of the goals of the 
present analysis, making use of a simple phenomenolog- 
ical parametrisation to estimate them. Moreover, the 
QCD DGLAP evolution can be introduced, which im- 
proves the data description in the large Q 2 kinematic re- 
ion accessible in the recent ZEUS DVCS measurements 
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This note is organized as follows. In the next section, 
we recall the main formulas for the color dipole formal- 
ism applied to DVCS. For the dipole cross section we have 
considered the saturation model [T^. which produces a 
unified and intuitive description of DIS [T^|. diffractive 
DIS vector meson production Drell-Yan polEU 
and DVCS 16] . In particular, the restriction to the trans- 
verse part of the photon wave function, due to the real 
final state photon in DVCS, enhances the contribution of 
larger dipole configurations and therefore the sensitivity 
to soft content and to the transition between hard/soft 
regimes. Such a feature provides a particularly relevant 
test of saturation models. Moreover, the approach in- 
cludes all twist resummation, in contrast with the lead- 
ing twist approximations. In the Sec. 3, we discuss the 
role played by the QCD evolution and skewedness in the 
high virtuality kinematic region. We also perform a sys- 
tematic analysis in order to investigate to what extent 
the distinct models improve the data description. These 
issues have implications in the correct determination of 
the t slope parameter £?, whose value has never been 
measured for DVCS. Finally, the last section summarizes 
our main results. 



II. DVCS CROSS SECTION IN DIPOLE 
PICTURE 

In the proton rest frame, the DVCS process can be 
seen as a succession in time of three factorisable sub- 
processes: i) the photon fluctuates in a quark-antiquark 
pair, ii) this color dipole interacts with the proton tar- 
get, iii) the quark pair annihilates in a real photon. 
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The usual kinematic variables are the -f*p c.m.s. en- 
ergy squared s = W 2 = (p + q) 2 , where p and q are the 
proton and the photon momenta respectively, the pho- 
ton virtuality squared Q 2 = —q 2 and the Bjorken scale 
x Bi =Q 2 /(W 2 +Q 2 ). 

The imaginary part of the DVCS amplitude at zero 
momentum transfer in the color dipole formalism is ex- 
pressed in the simple way [T^ . 

1 oo 

lmA(s,t = 0) = J dz J d 2 rH(z,r,Q 2 )a dlp (x,r 2 ) (1) 
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quark). The auxiliary variables e\ 2 = z(l — z) Q\ 
depend on the quark mass, to/. The K§ \ are the Mc- 
Donald functions and summation is taken over the quark 
flavors. 

Let us summarize the main features and expressions 
from the saturation model, which will be used here to 
estimate the DVCS cross section. A previous analysis 
compared to HI data can be found in Ref. [l^. The sat- 
uration model reproduces color transparency behavior, 
&dip ~ f 2 , for small dipoles, whereas it gives a constant 
behavior for large ones. This is rendered by a dipole cross 
section having an eikonal-like form, 



(To 



1 — exp 



) 2 

' sat 



(3) 

(*) = (f ) A GeV 2 , S = SB i (l + ^-) > (4) 



where the saturation scale Q sa t(x) (energy dependent) 
defines the onset of the saturation phenomenon and sets 
the interface between soft/hard domains. The param- 
eters were obtained from a fit to the HERA data pro- 
ducing g o = 23.03(29.12) mb, A = 0.288(0.277) and 
x = 3.04- 10~ 4 (0.41 • 10~ 4 ) for a 3-flavor (4-flavor) anal- 
ysis [T3. An additional parameter is the effective light 
quark mass, to/ = 0.14 GeV. For the 4-flavor analysis, 
the charm quark mass is considered to be m c = 1.5 GeV. 

The QCD evolution to the original saturation model 
was implemented recently [22J (BGBK), where the dipole 
cross section now depends on the gluon distribution in a 
Glaubcr-Gribov inspired way, 



g dip (x,r 2 ) = a Q 



1 — exp I — 



'' a s (p 2 )xG(x,iJ, 2 

3<70 



,(5) 
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H = i^r E e / { ^ + i 1 - z ) 2 ] £ i K ^ r ) 62 K ^ £2 r ) 
n f 

+ to 2 K (e ir )K (e 2 r)} , (2) 

where H(z, r, Q\ 2 ) = H* T {z, r, Q\ = Q 2 ) $ T (z, r, Q\ = 
0), with being the light cone photon wave function 
for transverse photons . Here, Q\ = Q is the virtuality of 
the incoming photon, whereas Q2 is the virtuality of the 
outgoing real photon. The longitudinal piece does not 
contribute at Q 2 = 0. The relative transverse separation 
of the pair (dipole) is labeled by r and z, (1 — z), are 
the longitudinal momentum fractions of the quark (anti- 
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FIG. 1: The DVCS cross section as a function of c.m.s. energy, 
W-yp. The curve is the result for the saturation model for fixed 
slope B = 4 GeV" 2 . 
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where the energy scale is defined as y 2 — C/r 2 
The parameters are determined from a fit to DIS data, 
with the following initial condition for LO DGLAP evo- 
lution, xG(x,fi 2 = lGeV 2 ) = A g x- X ° (1 - a;) 5 6 . The 
flavor number is taken to be equal to 3. The overall nor- 
malization (To = 23.03 mb is kept fixed (labeled fit 1 in 
Ref. )■ The DGLAP evolution improves the data de- 
scription in large Q 2 regime and brings the model close to 
the theoretical high energy non-linear QCD approaches. 

Having a suitable model for the dipole cross section, 
as in Eq. © or Eq. (J5J, we can use Eq. and then 
compute the final expression for the DVCS cross section 
as, 



1 [Zto.4(s,0)] 2 
a 7 p^ 1P ) = — — 

B 16 7T 



'1 



(6) 



where B is the t slope parameter (the behavior in |i| is 
supposed to obey a simple exponential parametrisation) . 

In our further calculations, the real part is included 
via the usual estimate p = tan(7rA/2), where A = A(Q 2 ) 
is the effective power of the imaginary part of the am- 
plitude. We have fitted it for 1 < Q 2 < 100 GeV 2 in 
the form A ofr (Q 2 ) = 0.2 + 0.01071n 2 (Q 2 /2.48). It can be 
verified that when it rises to ~ 0.3 at high virtualities 
the contribution of the real part can reach 20% of the 
total cross section. In the next section we compute the 
cross section above using the two versions for the satu- 
ration model and contrast them with the recent DVCS 
ZEUS data, which includes data points with larger Q 2 
values than the previous HI data. Moreover, we present 
a simple way to introduce skewedness effects into the cal- 
culation. 



III. 



RESULTS AND DISCUSSIONS 



In Fig. n is shown the result for the saturation model, 
Eq. ©, confronted to the experimental data on DVCS 
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FIG. 2: The DVCS cross section as a function of photon virtuality: (a) saturation model using B — 4 and 6.5 GeV -2 (solid and 
dashed curves) and Q 2 -dependent slope (dot-dashed curve - see text), (b) Effect of the BGBK model (includes QCD evolution) 
using B — 4 GeV -2 (solid and dot-dashed curves). 



of recent ZEUS measurements as a function of the c.m.s. 
energy, W lv (at fixed virtuality Q 2 — 9.6 GeV 2 ). The 
parameters of the 4-flavor fit have been used, producing 
good agreement with a fixed value for the slope, B = 4 
GeV~ 2 . 

In Fig. |21-a, we show the result of the saturation model 
for the behavior with Q 2 at fixed energy, W lv = 89 
GeV. In order to illustrate the sensitivity on the slope 
value, both values B = 4 GeV -2 (solid line) and B = 6.5 
GeV -2 (dot-dashed line) are shown 1 . Although the sta- 
tistical errors are large, it seems that for Q 2 >40 GeV 2 , 
the model underestimates the experimental data. This 
can indicate two things: (a) the slope diminishes as the 
virtuality increases or, (b) some additional effect appears 
at higher Q 2 . In order to investigate the first hypothe- 
sis, we compute the cross section using a Q 2 dependent 
slope, proposed in Ref. 23]. That is, B(Q 2 ) = B [I - 
0.151n(Q 2 /2)] GeV -2 which is based on the diffractive 
electroproduction of p. Such a slope dependence allows a 
good description of the Q 2 dependence of the cross sec- 
tion up to the highest measured values and gives a good 
normalisation for Bq = 5 GeV -2 . 

In order to investigate whether a QCD evolution im- 
proves the description, we show in Fig. |2}b the estimate 
using the BGBK dipole cross section, Eq. JSJ) as a func- 
tion of Q 2 using fixed slope values. There is an effect in 
the overall normalization and a slower decrease at large 
Q 2 in contrast with the model without QCD evolution re- 
producing well the ZEUS measurement for all Q 2 . This 
suggests that DGLAP evolution starts to be important 



[1] It is worth mentioning that a slope B = 6.5 GeV 2 (4-flavor) was 
able to describe correctly the HI experimental data for Q 2 < 40 
GeV 2 [T||. 



for the large Q 2 points measured by ZEUS. A comparison 
of the different Q 2 behavior independently of the normal- 
isation question is presented at the end of this section. 

Furthermore, we are motivated to investigate the im- 
portance of the skewedness effects in the DVCS process 
using the previous results. Here, we follow the approx- 
imation proposed in Ref. [24|. where the ratio of off- 
forward to forward parton distributions are obtained re- 
lying on simple arguments. The behavior of those ratios 
are given explicitly by p4| . 



Rq,g (Q 2 ) = 



2 2A+3 r (A + 



T(X + 3+p) 



(7) 



where p = for quarks and p = 1 for gluons, and where 
A is the exponent of the x~ x behavior of the input diag- 
onal parton distribution. It should be noticed that the 
skewed effect is much larger for singlet quarks than glu- 
ons. In the following, it will be assumed that the DVCS 
cross section is lead by a two gluon exchange. In our 
further computations, we use A = A(Q 2 ) as discussed in 
the previous section and the skewedness effect is given by 
multiplying the total cross section by the factor R 2 (Q 2 ). 
Once the effective power increases as a function of Q 2 , 
the skewedness effects could enhance the cross section by 
a factor two if values of A e // ~ 0.4 are reached at larger 
virtualities. In Fig. 0-a we show the result using the 
saturation model (4-flavor) and the skewedness correc- 
tion, Eq. {7|). The same analysis is shown for the BGBK 
model in Fig. |3J-b. The main effect is to increase the 
overall normalization of the cross section by about 40% 
and only slightly modify the large Q 2 behaviour. Again, 
this will be shown more clearly and independently of the 
normalisation at the end of this section. 

For completeness, we have investigated two additional 
versions of the implementation of the skewedness correc- 
tion factor. They are shown in Fig. 0] for fixed B = 4 
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FIG. 3: The results for the (a) saturation model with (full) and without (dot-dashed) skewedness effect and (b) BGBK model, 
with (dot-dashed) and without (full) skewedness effect for a fixed B = 4 GeV~ 2 . 
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FIG. 4: Comparison among different approximations for the 
skewedness correction (see text). 



GeV -2 . First, we have imposed skewedness correc- 
tion only for small dipoles by introducing R 9 {Xgbw = 
0.277) in the exponent of Eq. J3J| (dotted line). This 
is to prevent correction to the nonperturbative (large 
dipoles) piece of the dipole cross section. Further, we 
also test the rough approximation x = 0.41 XBj (dashed 
line), which comes from a simplified hypothesis a dip ~ 

Rg(X) (xbj) A . The conclusion is that these two differ- 
ent implementation of the skewedness correction do not 
make sensible changes w.r.t. the first skewedness correc- 
tion neither in normalisation nor in Q 2 dependence for 
the presently covered kinematic range and precision of 
the measurement. 

At this stage, some comments are probably needed. 
The estimate for skewedness taken into account above 
is an approximation as currently we have no accurate 
theoretical arguments how to compute it from first prin- 
ciples within the color dipole formalism. A consistent 



approach would be to compute the scattering amplitude 
in the non-forward case (the non-forward photon wave 
function has been recently obtained in Ref. 25]). In 
this case, the dipole cross section, <Jdip(%i, %2, t, A), 
depends on the light cone momenta x\ and xi carried 
by the exchanged gluons, respectively, and on the total 
transverse momentum transfer A ( additional informa- 
tion about the behavior on A is needed for the QCD 
Pomeron and proton impact factor) . The forward dipole 
cross section is recovered at x\ — X2 and A — 0. In the 
future, an experimental constraint for the nonforward 
dipole cross section should be feasible with increasing 
statistics on DVCS and exclusive (diffractive) vector 
meson production. 

To close this section, as the slope parameter B has 
never been measured for DVCS, we compare the differ- 
ent estimates presented in a systematic way separately 
for the effect on the Q 2 dependence and the effect on the 
overall normalisation. To compare the Q 2 dependences, 
we normalize all models to describe the ZEUS data point 
at the lowest Q 2 value, i.e. Q 2 = 7.5 GeV 2 . Further, we 
plot the ratio of each model to our baseline model SAT- 
MOD as a function of Q 2 . Such a procedure allows a Q 2 
dependence comparison independently of the normaliza- 
tion effect. These ratios are shown in Fig.[Sl-a, where the 
points (triangles-up) are the ratio of the ZEUS data to 
SAT-MOD including the error bars for the statistical (in- 
ner) and sum in quadrature of statistical and systematic 
(outer) uncertainties. 

On the other hand, to compare the effect on the nor- 
malisation we show the slope value needed to describe the 
lowest Q 2 value of the ZEUS data points B a = B(Q 2 — 
7.5 GeV 2 ). They are shown in Fig. 0-b. For complete- 
ness, we also present the measured slope values for vector 
meson production at that virtuality, both for p° and J/^ 
mesons as indications of typical values for respectively 
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FIG. 5: The ratio (a) MODELS/SAT-MOD as a function of Q 2 (models normalized to ZEUS data at Q 2 = 7.5 GeV 2 ), and (b) 
slope values of the models at Q 2 = 7.5 GeV 2 and slopes of heavy and light vector mesons. 



light and heavy mesons using the simple parametrisation: 



B = 0.60 
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(Q 2 



(8) 



where My is the meson mass. 

From these comparisons, we conclude that several 
models can account for the measured Q 2 dependence 
(SAT-MOD+B(Q 2 ), SAT-MOD+SKEW and BGBK, as 
well as combination of several of those effects) which 
are not distinguishable with the present experimental 
precision. The difference between the models is much 
more pronounced in the prediction of the cross section 
value, or in other terms, in the B value needed to de- 
scribe the integrated cross section over the available Q 2 
range. If the change in normalisation is small for the in- 
clusion of a Q 2 dependence in B, the effect is of the order 
of 12% for BGBK with respect to the basic SAT-MOD 
and of 40% for the skewedness effect (SKEW) and still 
larger when the different effect are combined (60% for 
BGBK+SKEW). 

In summary, these issues show clearly the importance 
of a measurement of the t slope parameter B. 



IV. SUMMARY 

It has been shown that the DVCS cross section at 
HERA can be described by the simple picture rendered 



by the color dipoles formalism. In particular, the satu- 
ration model does an excellent job in the current experi- 
mental kinematic domain. To achieve a good description 
of the data up the the highest Q 2 , the original saturation 
model can be supplemented by QCD evolution, an ad- 
ditional dependence of B on Q 2 and skewedness effects. 
These effects modify in a sensitive way the absolute cross 
section (10-60%). Measurement of the t slope param- 
eter B would already allow to discriminate among the 
different theoretical predictions with an amount of data 
comparable to the present ZEUS measurement. 



Acknowledgments 



M.V.T.M. thanks the CERN Theory Division, where 
part of this work was performed, for the hospitality and 
financial support. The work of L. Favart is supported 
by the FNRS of Belgium (convention IISN 4.4502.01) 
and M. Machado was partially financed by the Brazilian 
funding agency CNPq. 



[1] C. Adloff et al. [HI Collaboration], Phys. Lett. B 517, 
47 (2001). 

[2] L. Favart, Studies of DVCS and photoproduction of pho- 
tons at high t with the HI detector, In: International Eu- 



rophysics Conference on High Energy Physics (EPS03), 
July 17-23, Aachen (2003), hep-ex/0312013 . 
[3] S. Chekanov et al. [ZEUS Collaboration], Phys. Lett. B 
573, 46 (2003). 



[4] A. Freund, Eur. Phys. J. C 31, 203 (2003). 
[5] M. Diehl, Phys. Rept. 388, 41 (2003). 
[6] D. Miiller et ai, Fortschr. Phys. 42, 101 (1994). 
[7] X. Ji, Phys. Rev. D 55, 7114 (1997). 
[8] A. V. Radyushkin, Phys. Rev. D 56, 5524 (1997). 
[9] K. J. Golec-Biernat and A. D. Martin, Phys. Rev. D 59, 
014029 (1999). 

[10] L. Frankfurt et ai, Phys. Lett. B 418, 345 (1998), 
Erratum-ibid. 429, 414 (1998); A. Freund and V. Guzey, 
Phys. Lett. B 462 178 (1999). 

[11] J. C. Collins and A. Freund, Phys. Rev. D 59, 074009 
(1999). 

[12] A. Freund and M. McDermott, Phys. Rev. D 65, 074008 
(2002). 

[13] A. Freund and M. McDermott, Phys. Rev. D 65, 056012 
(2002); Phys. Rev. D 66, 079903 (2002); Phys. Rev. D 
65, 091901(2002); Eur. Phys. J. C 23, 651 (2002). 

[14] A. Donnachie and H. G. Dosch, Phys. Lett. B 502, 74 
(2001). 

[15] M. McDermott, R. Sandapen and G. Shaw, Eur. Phys. 
J. C 22, 655 (2002). 



6 



[16] L. Favart and M. V. Machado, Eur. Phys. J. C 29, 365 
(2003). 

[17] K. Golec-Biernat and M. Wusthoff, Phys. Rev. D 59, 
014017 (1999). 

[18] K. Golec-Biernat and M. Wusthoff, Phys. Rev. D 60 
114023 (1999). 

[19] A. C. Caldwell and M. S. Soares, Nucl. Phys. A 696, 125 
(2001). 

[20] M. A. Betemps, M. B. Gay Ducati and M. V. Machado, 

Phys. Rev. D 66, 014018 (2002). 
[21] M. A. Betemps, M. B. Ducati, M. V. Machado and 

J. Raufeisen, Phys. Rev. D 67, 114008 (2003). 
[22] J. Bartels, K. Golec-Biernat and H. Kowalski, Phys. Rev. 

D 66, 014001 (2002). 
[23] A. Freund, M. McDermott and M. Strikman, Phys. Rev. 

D 67, 036001 (2003). 
[24] A. G. Shuvaev, K. J. Golec-Biernat, A. D. Martin and 

M. G. Ryskin, Phys. Rev. D 60, 014015 (1999). 
[25] J. Bartels, K. Golec-Biernat and K. Peters, Acta Phys. 

Polon. B 34, 3051 (2003) 



